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A quench in an overdamped one dimensional cf> 4 
model is studied by analytical and numerical meth- 
ods. For an infinite system or a finite system with 
free boundary conditions, the density of kinks after 
the transition is proportional to the eighth root of 
the rate of the quench. For a system with periodic 
boundary conditions, it is proportional to the fourth 
root of the rate. The critical exponent predicted in 
Zurek scenario is put in question. 



Topological defects play a prominent role in many con- 
densed matter systems, see e.g. jjj for a review. It was 
also suggested that they were an important ingredient of 
the early universe S. Topological defects can be gener- 
ated in large numbers during a second order phase tran- 
sition. The dynamics of such a transition to a symmetry 
broken phase has been an object of much recent attention 
because of its importance in the cosmological context || 
and in condensed matter systems 

An early estimate of defect density after a quench was 
given by Kibble ||. It his theory, the speed of light 
is a dominant factor which determines the size of cor- 
related domains. An alternative scenario was put for- 
ward by Zurek j(| , who emphasized the importance of the 
nonequilibrium dynamics of the order parameter. Recent 
experiments point to the latter theory. The principal pre- 
diction of Ref. Q is the scaling law. He is able to predict 
by very general arguments an exponent in a power law 
dependence of the density of topological defects on a time 
of the quench in a variety of systems. The numerical ex- 
periment seems to confirm this prediction. 

It is normally assumed that close to the phase tran- 
sition the dynamics of the system is overdamped. We 
follow this line of argument to its logical limit and con- 
centrate from the very beginning on a purely diffusive 
dynamics. The model we study is the <^> 4 theory in one 
spatial dimension. Although there is no proper phase 
transition in such a one dimensional model there is a 
long range order in its symmetry broken phase. The real 
field 4> falls into one of the two minima of the Ginzburg- 
Landau potential at0 = +lor0 = — 1. The domain 
walls separating these two different vacua are the topo- 
logical solitons or the kinks. We develop a novel method 
of how to predict the density of kinks after a quench. It 



turns out that the density is proportional to the eighth 
root of the quench rate in contradiction with Zurek's pre- 
diction of the fourth root. We point to an implicit as- 
sumption in Zurek scenario, which lays at the bottom of 
this discrepancy. 

Let us consider a one dimensional overdamped 4> 4 
model 

<p(t, x) = (j>"(t, x) - 2a(t)(j)(t, x) - 2<p 3 (t, x) + rj(t, x) , 

(1) 

where ' = dt, ' = d x . 4>(t,x) is a real order parameter. 
r)(t, x) is a Gaussian noise of temperature T with corre- 
lations 

< rj{t,x) >= , 

<r 1 (t 1 ,x 1 )r 1 (t 2 ,x 2 )>=2T6(t 1 -t 2 )6(x 1 -x 2 ) . (2) 

The coefficient a(t) is time dependent. We consider a 
linear quench 



A , if t < 

a(t) = i A(l - i) , if < t < r^±i 
-1 ,if r^<t 



(3) 



Before the quench, for t < 0, the system is in a symmetric 
phase (A > 0), during the quench, at t = r, it undergoes 
a transition from the symmetric phase a{t < r) > to a 
broken symmetry phase a(t > r) < 0. Finally it settles 
down at a(t) = — 1. 

As long as the system is still in the symmetric phase 
(t < t) , the field <j) can be regarded as a small fluctuation 
around its symmetric ground state <p(t, x) = 0. It is 
justified to neglect the cubic term on the RHS of Eq. ([!]). 
The field can be written as a Fourier transform, 
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/+oo 
dk 4>{t,k)e lkx . 
-oo 

The Fourier transform of linearized Eq. (|IJ) is 

<t>{t, k) = ~[2a(t) + k 2 )4>(t, k) + fj(t, k) . 
The Fourier transformed noises have correlations 
<fj{t,k) >=0 , 

T 

< rj*(ti, kx)rj(t 2 , k 2 ) >= -<% - t 2 )8(k 1 - k 2 ) . 

IT 

The formal solution of Eq. (||) at the time r is 



(4) 



(•5) 



(6) 



1 



<j>(T,k)= f dtx exp{- I dt 2 [2a(t 2 ) + k 2 ]}fj(t,k) . 

(7) 

This solution, the correlations (^|) and the explicit form 
of a(t), see Eq.(||), give the correlation function at t = r 

< 4>*{t 1 M)kt2,k 2 ) >= -S(h-k 2 ) G(t, A, k) 

7T 



G(r,A,k) 



-2k 2 r 



AA + 2k 2 



'^f^-^A + k^-Erf^-k 2 )] . (8) 

The linearized approximation can not be used in the bro- 
ken symmetry phase. Nevertheless we would like to use 
the correlations at t = t to predict the density of topo- 
logical defects after the quench. We have to make con- 
tact between the formalism developed for the symmetric 
phase and the one which is relevant to the symmetry 
broken phase. 

a(t) = — 1 in the symmetry broken phase. If the noise 
is neglected, then Eq.(|l|) has a static kink solution 



cj>(t, x) = tanh(a; — z) = K(x — z) 



(9) 



with an arbitrary kink position z. The field which is 
going to relax to this equilibrium configuration (|^) can 
be expressed as a sum of the kink solution and of eigen- 
functions u a of the fluctuation operator around the kink 
numbered by a discrete/continuous index a, 

(j>{t,x) = K(x - z) +^2<Z> a (t)u a (x - z) . (10) 



The eigenstates have positive eigenvalues and their am- 
plitudes Q a (t) decay exponentially with time. The "kink 
transform" of the field (Coh 



I(t,x) = 



dyK'(y)(f>(t,x + y) 



(11) 



vanishes for x = z because K'(x) is orthogonal to K(x) 
and to u Q (x)'s. The zero of the kink transform of the 
initial field is the equilibrium position of the kink. The 
same transform ( p"l] ) can be used to predict positions of 
antikinks. This method has been introduced and exten- 
sively tested numerically in g]. To predict the positions 
of kinks after the quench one must find the kink trans- 
form 7(r, x) of the initial field 0(r, x) and locate the zeros 
of 7(t, x). 

We do not need any detailed knowledge about the ac- 
tual positions of kinks. All we need to know is their 
average density. The density of kinks n can be expressed 
by the density of zeros of the kink transform 7(t, x) 



<N> 
2L 



= h f L dx <IA^)Wt,z)]> 



(12) 



L is the length of the considered interval. Due to trans- 
lational invariance 



n =< \I'(t,0)\ <5[/(r,0)] > 



(13) 



Thanks to the same symmetry < 7(r, x)I'(t, x) >= for 
any x, so that 



n=< \I'(t,Q)\ ><5[I(t,Q)] > 



(14) 



The identity |7'| = I' Sign(I') and the Fourier trans- 
forms of the (5 and Sign functions help to derive 



7T / < 7' 2 (t,0) > 

71 - 2 V <72(r,0) >■ 



(15) 



The higher is the ratio of the variation of 7 to its average 
magnitude, the higher is the density of kinks. The higher 
is this ratio, the more likely is the function 7 to cross zero. 

One can express the kink transform with the Fourier 
transform of </>(r, x) (0) as 



dk U(k) 4>(T,k) , 

■OC 

dk k U{k) 4>(t, k) 



7(r,0) 
7'(r,0)=i 
U(k) = - 



■nk 



sinh(7rfc/2) 



(16) 



U{k) is a Fourier transform of K'{x). The 7 — 7 correla- 
tions in Eq.(|l5|) can be expressed by the correlations (g), 
so that finally we obtain 



7T Jdk k 2 U 2 (k)G(T,A, k) 
n ~ 2 V JdkU 2 (k)G{T,A,k) 



(17) 



The kink size provides an ultraviolet cut-off through the 
function U 2 (k). Without any UV cut-off the expression 
( |l7| ) would be infinite for any A and r. In the symmetric 
phase the noise driven order parameter <p(t, x) has an 
infinite number of zeros. 

To extract the scaling law, we must find the asymp- 
tote of n for large r. The asymptote of the correlation 
function G is 



G(t , A, k) 



T e 2 Ar+^ / dxe - 

2^ I/*:* 



(18) 



We can make furthe r ap proximation under the integral 
in the numerator of (|17|) 



r 



dk k 2 U 2 {k) G(r,A,k) » 

/ + OO A 
dk [/2 (fc) = ^ e 2Ar 
-oo 3 



(19) 



We introduce p — k i^) 1 ^ 4 in the denominator integral 
to obtain the asymptote 



2 



p +00 



dk U 2 {k)G(r,A, k) 
+00 



2A^ 



dp U 2 



+ OO 



+00 



(20) 



After the numerical coefficients are worked out one ob- 
tains 



1.16 



(21) 



The density of kinks is proportional to the eighth root of 
the transition rate. The asymptote is achieved for r >> 
2A. The critical exponent does not depend on the actual 
choice of the UV cut-off provided by the function U 2 (k). 
In a lattice theory regularized by its lattice constant, 
one would get the 1/8 critical exponent for the density 
of zeros of the order parameter <p(r, x) itself. 

The last prediction must be contrasted with a corre- 
sponding result for a model with periodic boundary con- 
ditions. The periodicity <f>(t, L) — <p(t, 0) implies that 
momentum is quantized k m — (^-)m with integer m. 
The integrals over k in (|l7| ) are replaced by summations 
over m and all appearances of k are replaced by k m . The 
denominator sum is dominated for large r by the m = 
mode contribution, 



E U 2 (k m )G( T ,A,k m ) 

a— — 00 
+ 00 I 

E U 2 (k m )J^-e 2A ^ 



+ OC 



dx e 



«-|-oo 



'— U 2 (0) 



(22) 



The numerator sum is approximated in a similar way as 
the numerator integral for the infinite system, 



E 



+00 

k 2 m U 2 (k m ) G(r, A, k m ) 

m— — 00 
+00 

T 2A-H 



E k 2 m U 2 (k m ) x /—e 



m— — 00 



+00 



dx e 



T 1,2 



E ^ 



(23) 



Putting these two asymptotes together we get 



(24) 



If the periodic boundary conditions are imposed, the den- 
sity of kinks is proportional to the fourth root of the 



transition rate. In fact this conclusion is the same for 
any fixed boundary conditions which quantize k. The 



(2t) 



» 1. 



asymptote is achieved for 



. 2/1 '"1 V 2A 

We expect for a large but finite system, that the density 
scales like the eighth root of the quench rate for moder- 
ately large r but for very large r it scales like the fourth 
root. Transition between these two regimes takes place 
at a r which grows like the fourth power of the system 
size L. The thermodynamic limit and the large quench 
time limit do not commute. Let us note in passing that a 
similar dependence of critical exponent on boundary con- 
ditions takes place in the diffusion of overdamped kinks 
in the sine-Gordon chain ||. 

The Zurek scenario as described in || does not refer 
to boundary conditions. Let us analyse our derivation 
from the point of view of Ref . || . The relaxation time of 
the system depends on time like t r = 2 a(t) • ^ ne system 
ceases to follow the changes in a(t) at the instant t, when 
its relaxation time becomes greater than the time left to 
the transition at t = r, (t — t) — - \\ . For the linear 

v ' 2a[t) 

quench (^) this takes place when (r — t) = \l At this 
instant of time the correlation length is £ = (■^) 1 ^ 4 - The 
density of kinks after the transition is inversely propor- 
tional to this "frozen" correlation length, so n ~ (^) 1//4 . 

The weak point of this argument is that t r = 2 a(t) 
is the relaxation time but only for the k = mode. 
Modes with nonzero momenta k > have relaxation 
times t r (k) = 2 „.(£\+k i ^ ^"(0) and they are frozen later 

As we can see in Eq.(Ea), the 



2a(t)+k' 2 

than the k — mode 
density of nucleated kinks is related to the ratio of the 
magnitude of short wavelength modes to the magnitude 
of long wavelength modes. Contrary to Zurek scenario, 
which implicitly assumes that all the modes are frozen at 
the same instant of time, the k > modes can still grow 
for a while even after the k = mode was already frozen. 
Thus the density of kinks after quench is greater than 
that given by Zurek estimate. It is a matter of detailed 
calculations presented above to find out that the actual 
critical exponent is 1/8 instead of 1/4. 

It is the point to discuss the numerical experiment of 
Laguna and Zurek [Q. Those authors discuss the model 
(|l|) with (f> added on the LHS. Their numerical experi- 
ment on a lattice with periodic boundary conditions gives 
a crititical exponent of 1/4. We could naively argue that 
this is due to the periodic boundary conditions and in 
agreement with our theory. This is, however, not the 
case. The lattice is large enough to render the bound- 
ary conditions irrelevant. Hence according to our the- 
ory one would expect the 1/8 scaling. But the point 
is that neither Zurek estimate nor our calculations ap- 
ply to this system. The damping coefficient 77 = 1 used 
in actual simulations (in distinction to the theoreti- 
cal discussion |Q) is in the underdamped range. At the 
beginning of the quench all the field oscillators with dif- 
ferent momenta k are underdamped. The k = would 
be overdamped, if 7/ > \/2 in units of Ref. 0. The short 



3 



wavelength modes would require even higher r\. At the 
beginning of the quench the second order time derivative 
dominates over the first order derivative. Without any 
reliable tool at hand we can only make a very superficial 
remark. Replacement of 4> by <j) should give the replace- 
ment of t by r 2 in any scaling behavior. This replacement 
in our formula leads to the critical exponent 1/4 in 
agreement with the numerics of Ref . . 

We performed a numerical simulation of the quench 
(||) in the model ([!]). The total extension of our lattice 
was 20 units, the lattice constant 0.1 and the temper- 
ature T = 10~ 8 . With these parameters the periodic 
lattice system was far in the thermodynamic regime, the 
asymptote with the 1/4 scaling is expected to be achieved 
for r >> 10 4 . The parameter A in Eq.([|) was chosen as 
A = 100. We have measured the density of zeros of the 
function 0(r, x) for r = 0.001 . . . 3.2 units of time. The 
simulation was repeated 20 times for each r. The fig- 
ure shows the double logarithmic plot of the density of 
zeros n z as a function of the quench time. The slope is 



of 



0.110 — 0.005, which is consistent with our prediction 
i = -0.125. 



Conclusion. 

The critical exponent observed in an experiment on 
a macroscopic system in its overdamped regime should 
be 1/8. One can observe the 1/4 scaling in a numerical 
"experiment" on a small lattice. 
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